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Modification of the renormalization-group approach, invoking Stratonovich transformation at each 
step, is proposed to describe phase transitions in 3D Ising-class systems. The proposed method is 
closely related to the mean-field approximation. The low-order scheme works well for a wide thermal 
range, is consistent with a scaling hypothesis and predicts very reasonable values of critical indices. 
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It is well understood since 1960's that an essential in- 
crease of fluctuations accompanies a second-order phase 
transition. Exact results for 2D Ising model 0] showed 
that Landau phenomenology Q was too crude to describe 
the vicinity of the transition point. Ginzburg demon- 
strated that it is correlations what breaks Landau ap- 
proach down ||. He established the parameter which 
determines the thermal range for the so-called fluctuation 
region. According to the scaling hypothesis H , the prop- 
erties of a system in the fluctuation region are determined 
by the single quantity r c , being the correlation length for 
fluctuations. This parameter shows a power-law depen- 
dence on the external field h and on the thermal interval t 
from the transition: r c cx , r c cx t~ v . Values of and 
v are determined by the universal characteristics, like the 
dimensionality D, the number of components of the order 
parameter, and the type of the interaction in the system 
(long-range or short-range). Other critical indices can be 
expressed via fj, and v. For instance for the scalar order 
parameter x one obtains x cx , f3 = vD — v / ' /.i. 

The general framework to obtain critical indices is the 
renormalization-group (RG) approach § |, §. The RG 
transformation consists of the two stages. First, one in- 
tegrates out all Xk with wave vectors k > fco/s, where ko 
is a cut-off and s > 1 is a parameter. The main assump- 
tion of the method is that the potential energy of thus 
renormalized system can be truncated to certain simpli- 
fied form, for example to the form of <fi 4 model. The 
second step is a scale transformation. It is aimed to re- 
store the cut-off for k back to fco and preserve the original 
dispersion law for the long-wavelength modes. 

The RG transformation should be applied recursively; 
its unstable fixed point corresponds to a phase transition. 
Lyapunov exponents for the vicinity of the fixed point 
determine the critical indices. The scaling hypothesis 
requires exactly two of the Lyapunov exponents to be 
positive, they correspond to the two independent critical 
indices. 

With respect to the value of s, two modifications of 
the RG scheme can be distinguished. Wilson's e expan- 
sion H ^| deals with s> 1. The method is formally valid 
near the upper critical dimension D up (for the Ising-class 
classical systems D up — 4). The "most divergent" dia- 



gram series can be pointed out in this case. This gives 
the asymptotic expansion in terms of e = D up — D. In 
the continuous version of the scheme (so-called exact RG 
approach) |Q, s is infinitesimally close to unity. Certain 
decoupling (truncation) should be introduced in this case 
to solve the integro-differential equation for RG flow. An 
accuracy is then determined by a form of the decoupling 
used. 

All the approaches mentioned above are very well- 
developed. They allow to calculate critical indices and 
universal combinations of amplitudes with a high accu- 
racy. Progress is also achieved in the calculation of non- 
universal quantities; particularly a crossover between 
Gaussian (Landau-like) and critical behaviour is analyzed 
B. However RG is applicable only for long- wavelength 
excitations since they can be described as a fluctuat- 
ing field. Therefore, in a practical calculation, short- 
wavelength excitations should be first integrated out by 
means of some perturbation expansion, and then RG 
should be applied. The short- and the long-wavelength 
excitations should be separated ad hoc. Such a two-step 
procedure is not very elegant. At least from the method- 
ological point of view it is desired to develop a scheme 
which gives all the information in a uniform manner. 

It is interesting to note that calculations are much sim- 
pler outside the critical region. In fact, even the mean- 
field approximation (MFA) is usually sufficient for the 
qualitative analysis of this region. In MFA, the non-local 
interaction in the system is replaced by an effect of a 
self-consistent average field. The MFA predictions for 
the transition temperature and thermal behaviour of the 
order parameter in 3D case are quite reasonable |I. On 
the other hand MFA gives Landau set of indices and is 
not very accurate. Note that Landau approach itself has 
a mean-field nature. 

In this paper, we present a simple generalization of 
MFA, performing well for the 3D Ising-class systems both 
nearby and far from the transition. Particularly the 
critical behaviour is reproduced correctly: the values of 
indices are predicted with a few-percent accuracy. Es- 
sentially, the present approach consists in the alternat- 
ing application of the Stratonovich transformation and 
the most primitive version of the renormalization-group 
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transformation at a finite value of s. 

We consider a classical scalar field with an anharmonic 
local potential and a short-range harmonic interaction at 
temperature T. The partition function of this system is 

Z = J [Dx]e- V ' T - (1) 

V = [u(x r + x)dr + J2^^- 
J k 2 

Here U(x-\-x) is an even on-site potential, x k is a Fourier 
transform of the displacement x r from the average posi- 
tion x; x ^ below the transition temperature, and the 
average of x r equals zero by definition. To make the 
definitions unambiguous, we put additional conditions 
U(x) = 0; f2fc=o = 0. The wave vector k does not exceed 
the unititary cut-off: k 2 < 1. The case of a short-range 
interaction corresponds to the asymptotic behaviour 

n k =u;k 2 ,k^0. (2) 

The simplest way to integrate out the modes is to skip 



Xk with k > 1/s. This corresponds to the zero-th approx- 
imation in the nonlinearity of U. The subsequent scaling 
results in the RG transformation k — > sk, U — > s 3 U, x — > 
x/y/s. The defect of thus applied RG method comes 
from the the increase of correlations (non-linear part of 
U grows). Therefore, after several iterative transforma- 
tions, the RG procedure is not valid any more. Although 
the formal analysis of the fixed point U(x) — can be 
done, it wrongly gives Landau-like critical indices. 

To avoid the problem higher-order approximations and 
the truncation of the potential are used in known schemes 
described above (a very clear explanation is given in |fTof ). 
Here we use another procedure: the zero-th approxima- 
tion is used, but the change of variables is made at each 
RG step. The change of variables reduces correlations 
and is aimed to compensate its increase at RG transfor- 
mation. 

At the first stage, we integrate out (omit) all the modes 
with £l k > fj, where 0, = J fl k d 3 k is an average of il k . 
For the rest of the modes, the Stratonovich identity is 
utilized: 



exp 



(n k -Ci)\x k \ 2 

T 



1T& 



* j exp (-1 hme^kfu) 



k 



dfk 



(3) 



Here the integration over the complex variable f k — $l{fk) + z3(/fc) denotes the integration over the complex plane: 

After (||) is substituted into (|l|), expressions J2 k ^- e ( x kfk) = / x r f r dr and J2k^\ Xk \ 2 = I ^ x r can be collected in 
the exponent. Therefore 
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where f r is the inverse Fourier transform of f k ■ 

Now, we can integrate over x r . Let us introduce the function F(f) accordingly to the equation: 
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where Fq is a constant, defined from the condition F(x) = 
0. 

With thus defined F, the partition function takes the 
form 



Z (x 
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(6) 



W = f F{f r + x)dr + J2 



nn k \f k \ 2 
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Equation (||) is formally very similar to Eq. ([!]). Func- 
tion F is even; the dispersion at k — > coincides that 



of Eq. ([!]). One can also guess from Eq.(||) that as the 
average of x r equals zero, the average of f r also vanishes. 
Indeed, if this average would not be zero, f k= a would 
fluctuate around a macroscopically large average value, 
and (0) would not be fulfilled. 



We argue here that variables f k are much less corre- 
lated, than x k . Let us consider the parabolic approxima- 
tion for F: F(x + f) w ^/ + ^gf 1 / 2 . Variables f k 
are uncorrelated in this approximation. Obviously < / > 



3 



vanishes only if 



dF(x) _ 
df 



= 0, that gives 



fix _ J x exp (-T- 1 (U(x + x) + \Ftx 2 )) dx 
~T ~ J exp (-T- 1 (U(x + x) + ^Clx 2 )) dx ' 



(7) 



This is exactly the mean-field equation of state for the 
system (Q). So, even the complete neglect of correlations 
in /fc still gives something reasonable. On the other hand, 
the parabolic approximation for U directly in Eq. (|l|) gives 
non-sense: since U(x) does not depend on temperature, 
the system does not show a transition at all in this ap- 
proximation. Such a comparison convinces us that cor- 
relations in fk are less essential, than in Xk- As it was 
pointed above we expect that this can compensate an 
increase of correlations in RG transformation. 

Now, we finish the RG transformation, operating with 
Eq.(^|). All fk with k > 1/s are to be integrated out 
(omitted in our approximation) . After the scale transfor- 
mation fc — > sk, / — > x/^/s is done, the potential energy 
of the system takes the form 



where 
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FIG. 1: Upper panel: the fixed point of the transformation 
( p"o| ) at s = 2.5. Lower panel: critical index v vs. scaling 
factor s. Note a very weak dependence of v(s). 



C/ (1) (x + 5 (1) ) = s 3 F(x/^ + x), x (1) =x^. (10) 

Such a transformation should be applied recursively, 
like in the standard RG approach. The analysis can be 
easily performed numerically. 

The sequence fi/c, ^!?\ ■■•is independent of U and 
T. At every step a new fl^ are to be calculated. After 
several transformations and converge to certain 
limit. Calculation from formula (m gives 
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is implicitly defined from the equation Q} 



oo) 



u^d*k. 



Now, let us consider the properties of the sequence 
U, U^ll/W,... at x = 0. Transformation -> U {n+1) 
has a non-trivial unstable fixed point Uf(x). This func- 
tion is shown in the upper panel of Fig.l. Consider small 
deviations from this fixed point: U^ n \x) — Uj(x) + 
u n (x). The linearization of the transformation (nil) gives 

u n+ i{y) = u° + s 3 u n (x)A(x,y)dx; (12) 



exp 
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where u° stands to fulfill the condition u n+1 (0) — 0. The 
analysis of this operator shows that among its eigenvalues 
only two exceed unity. Denote them s Ac and s A ° ; they 
correspond to an even and odd eigenvector, respectively. 
Thus the theory is consistent with the scaling hypothesis, 
as it was mentioned in an introductory part. 

Exponents A e and A D are the inverse of the critical 
indices v and /i, respectively I], 10 . In the ideal case, 
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critical indices should be independent of the particular 
U(x) and flk, as well as of the value of the factor s. 
Indeed, after a proper re-scaling of the units, the values 
of T and drop out from Eq.(10). Therefore in our 
model v and /i can depend only on s. Further, it turns 
out that n = 2/5 at any s, as it should be for the local- 
potential approximation Q. What is about v, it depends 
on s very weakly. The numerical result for v(s) is shown 
in the lower panel of Fig.l. 

On the other hand it is important to point out that the 
crude universal properties of the model enter the trans- 
formation (|To| ) essentially. The behaviour flk^o k 2 
corresponds to a short-range range interaction. Such be- 
haviour and the dimensionality of the model determine 
the scale transformation, i.e. powers of s in Eq.©. It 
is also directly reflected in the formula that the order 
parameter is a scalar. 

For the more detailed analysis it is reasonable to con- 
sider an extreme value of v{s), because in this case small 
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TABLE I: Values of the critical indices obtained from the 
Lyapunov exponents for the transformation ([l2]) at s — 2.5. 
Results of the Landau theory and known values || for 3D 
Ising model are given for the comparison. 

Index Estimation from Eq.(^) 3D Ising Landau theory 
v 0.652 0.63 1/2 
(i 2/5 0.403 1/3 
P 0.326 0.327 1/2 



variations of s do not affect v. This gives s rs 2.5, v 
0.652. The obtained values of indices are collected in Ta- 
ble 1. An agreement with the numerical results for the 3D 
Ising model is even better than one could expect. Since 
the scaling hypothesis holds in the model, the values of 
all other critical indices can be determined and should 
have a small error bar. As an example, the calculated 
value of f3 is presented in Table 1. 

The advantage of the present approach is that it allows 
to estimate not only critical indices, but also the (non- 
universal) thermodynamic quantities in a wide thermal 
range. Here we calculate the dependence of the order 
parameter x on temperature for the 3D Ising model. The 
model is defined as a discrete cubic lattice with U (x) — 
5(x— 1) + S(x + 1) and the nearest-neighbors interaction 
Ofc = 3 — cos nk x — cos nk y — cos irk z . The initial Brillouin 
zone is a cube —1 < k Xt y tZ < 1; at the first step the value 
ft = 3 is used and all the modes with |fc| > 1/s are 
integrated. The procedure for further steps is exactly 
as described above. To simplify the calculation, we put 
Qy 1 ) — fl(°°> at these steps; this simplification almost 
does not change the result obtained. To determine the 
value of x we use the condition 



d . . 
-j^U n {x n ) = Q,n 



(13) 



(compare with the derivation of formula (0)). Numeri- 
cally, it is enough to calculate U n up to n m 8. The result 
for x 2 (T) is presented in Fig. 2. The mean-field and ex- 
perimental (Monte-Carlo) data are sketched in the same 
figure for comparison. Far from the transition is al- 
most parabolic, therefore the theory passes into MFA. 
The critical behaviour occurs at the transition point. 
Note that Ising model is very non-linear, but the scheme 
still performs good. 

We conclude that the presented low-order scheme de- 
scribes correctly all the physics of the second-order phase 
transitions in 3D systems with a short-range interaction 
and a scalar order parameter. It gives a unified descrip- 
tion both of the fluctuation region and of the system far 



from the transition point. The method does not pretend 
to be very accurate; its advantage consists in simplicity. 
In principle, one can improve the accuracy by a higher- 
order procedure for the integration in RG transformation. 
An applicability of the method for another universality 
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FIG. 2: The dependence of the square of the order parameter 
on the temperature in the 3D Ising model. Solid curve: the 
present approach with a — 2.5; the dot curve: mean-field 
calculation; the dash curve: numerical results. Thin ticks 
show the tangents to x 2 (T) in the transition point. 



classes is the subject of further studies. 
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